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l3 preliminary data of two-particle Bose-Einstein correlations are reported for hadronic 
Z° decays in c + c~ annihilation at LEP. The invariant relative momentum Q is identified 
as the eigenvariable of the measured correlation function. Significant anti-correlations are 
observed in the Bose-Einstein correlation function in a broad region of 0.5 — 1.6 GeV with 
a minimum at Q s» 0.8 GeV. Absence of Bose-Einstein correlations is demonstrated in the 
region above Q>1.6 GeV. The effective source size is found to decrease with increasing value 
of the transverse mass of the pair, similarly to hadron-hadron and heavy ion reactions. These 
feautes and our data are described well by the non-thermal r-model, which is based on strong 
space-time momentum-correlations. 
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Introduction: Boson interferometry provides 
a powerful femtoscopic tool for the investigation 
of the space-time structure of particle produc- 
tion processes on the femtometer lengthscales. 
Bose-Einstein correlations (BEC) of two identi- 
cal bosons reflect both geometrical and dynamical 
properties of the particle radiating source [l|-0] • 

In e + e~ annihilation BEC have been ob- 
served 0] to be maximal when the invari- 
ant momentum difference of the bosons, Q = 
a/— {pi — P2) 2 , is small, even when one of the 
relative momentum components is large. This 
is not the case in hadron-hadron interactions [j| 
or in heavy-ion interactions [lol j. where BEC are 
found not to depend simply on Q, but to decrease 
even if Q is small but any of the relative momen- 
tum components is large. The size (radius) of the 
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source in heavy-ion collisions has been found to 
decrease with increasing transverse momentum, 
Pt, or transverse mass, m t = y m 2 + p 2 , of the 
bosons, for a recent on data oriented review see 
ref. [f|. A similar effect has been seen in p+p col- 
lisions 11], as well as in e + e _ annihilation [r|. 
Such a behavior that can be described by hydro- 
dynamical models of the source, 0], however, 
a simple Q dependence of Bose-Einstein correla- 
tions in e + e~ collisions can not be exp lained in 
hydrodynamical or thermal models [131 ] . 

Event and track selection: The data used in 
the present analysis were collected by the L3 de- 
tector at LEP at an e + e~center-of-mass energy 
of y/s ~ 91.2 GeV. In total about 0.8 million 
events with an average number of about 12 well- 
measured charged tracks are selected. This re- 
sults in approximately 36 million like-sign pairs 
of well-measured charged tracks. Events are clas- 
sified as two- or three-jet events on the basis of 
the Durham jet algorithm with a jet resolution 
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parameter y cut = 0.006, yielding about 0.5 mil- 
lion two-jet and 0.3 million three-jet events. In 
the present study we report only about the results 
for two-jet events, more detailed three-jet results 
will be presented elsewhere. 

Bose-Einstein correlation function: The 
two-particle correlation function of two particles 
with four-momenta p± and p2 is given by the ra- 
tio of the two-particle number density, P2(pi,P2), 
to the product of the two single-particle num- 
ber densities, pi(pi) Pi(p2) ■ Since we are here 
interested only in the correlation R2 due to 
Bose-Einstein interference, the product of single- 
particle densities is replaced by po(pi,P2), the 
two-particle density that would occur in the ab- 
sence of Bose-Einstein correlations: 



R2(Pl,P2) 



P2(Pl,P2) 

Po{pi,P2) 



(1) 



An event mixing technique is used to construct 
P0j whereby all tracks of each data event are re- 
placed by tracks from different events having sim- 
ilar multiplicity to the original event. 

po is corrected for detector acceptance and 
efficiency in the same way as p2- The mixing 
technique removes all correlations, for example, 
resonances and energy- momentum conservation, 
not just Bose-Einstein correlations. Hence, po 
is also corrected for this by a multiplicative fac- 
tor which is the ratio of the densities of events 
to mixed events found using events generated by 
Jetset without BEC simulation. Thus R2 
is measured by 

R2 = {R2 data-^2 gen) / (-^2 det-^2 gcn-noBE ) j \^) 

where data, gen, det, gen-noBE refer, respec- 
tively, to the data sample, a generator-level 
Monte Carlo sample, the same Monte Carlo sam- 
ple passed through detector simulation and sub- 
jected to the same selection procedure as the 
data, and a generator-level sample of a Monte 
Carlo without BEC simulation. 

The invariant relative momentum Q is 
eigenvariable of the correlation function in 
e + e~ annihilation. In e + e~ annihilation at lower 
energy [1] it has been observed that Q is the 
appropriate (eigen)variable of the Bose-Einstein 
correlation function, which implies an approxi- 
mate spherical symmetry of particle emission in 



the rest frame of the pair. A priori, one does 
not expect the hadron source to be so spheri- 
cally symmetric in jet fragmentation. Recent in- 
vestigations have, in fact, found an el ong ation of 
the source along the jet axis [12, Il5l4l8]. While 



this effect is well established, the elongation is 
actually only about 20%, which suggests that a 
parametrization in terms of the single variable Q, 
may be a good approximation. 

We have checked on l3 data, if indeed Q is 
an eigenvariable of the BEC or not, and confirmed 
[l9| that this is indeed the case, both for all and 
for two-jet events: We observe that R2 does not 
decrease when both q 2 = (p\ — P2) 2 and q^ = 
(E\ — E2) 2 are large while Q 2 = q 2 — q$ is small, 
but is maximal for Q 2 = q 2 ' — q^ = 0, independent 
of the individual values of q and go- Furthermore, 
two-dimensional fits with the parametrization 

R 2 (q, qo) = 1 + A exp ((rq) 2 - (r q ) 2 ) (3) 

find r and ro to be equal. (Where we also note 
parameter A, the intercept parameter of the corre- 
lation function.) The similar conclusion is found 
in a different decomposition: Q 2 = Q 2 + Q\ B , 

where Q 2 = (pn — Pt2) 2 is the component trans- 
verse to the thrust axis and Q\ B = (p\\ — pn) 2 — 
(E\ — E2) 2 combines the longitudinal momentum 
and energy differences. Again, R2 is maximal 
along the line Q = 0. This is a non-trivial re- 
sult. For a hydrodynamical type of source, on 
the contrary, BEC decrease when any of the rel- 
ative momentum components is large @, H|. 

In the region of < Q < 0.5 GeV, we ob- 
serve as usual, a positive correlation, due to the 
Bose-Einstein symmetrization effect of like-sign 
charged identical boson pairs. In the region 0.5 
< Q < 1.6 GeV, the measured l3 correlation 
function R2 decreases below unity [3l[, which is 
indicative of an anti-correlation. This is clearly 
seen in Fig.[T]by comparing the data in this region 
to an extrapolation of a linear fit, 7(1 + eQ) that 
is fitted to our data in the region Q > 1.6 GeV, 
where 7 as an absolute normalization constant 
and e is a measure of long-range, residual non- 
Bose-Einstein correlations in our measurement. 
The extrapolation to the low values of the in- 
variant relative momentum Q is indicated by the 
dashed line on Fig.[TJ Correlation functions with 
1 + positive definit forms are by definition un- 
able to describe this dip in R2. This is the pri- 
mary reason for the failure of several Q dependent 
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FIG. 1: l3 data on Bose-Einstein correlations of 
charged particles are compared to a r model fit. The 
fit describes the peak structure at low Q < 0.5 GeV, as 
well as the region of anti-correlations in 0.5 < Q < 1.5 
GeV around Q w 1 GeV with a x 2 /NDF = 90.1/95. 
Note that the background distribution in the Q > 2 
GeV is within errors flat, all long-range correlations 
have been removed. A indicates theory - data. 



parameterizations. We note that this dip is less 
apparent, if one only plots (and fits) R 2 for Q < 2 
GeV as has usually been done in the past. 

Many parametrizations discussed earlier, for 
example Gaussian or Levy source distributions as 
well as Edgeworth expansion, have been shown 
before to be insufficient to describe the BEC (l9l - 
211 ] . These parameterizations assume a static 
source: the parameter R, representing the size 
of the source as seen in the rest frame of the pion 
pair, is a time independent constant. It has, how- 
ever, been observed that R depends on the trans- 
verse mass, mt = \J m 2 + p\ = y/ ' E 2 — p 2 , of the 
pions [12] • It has been shown [2^, [23|] that this de- 



pendence can be understood if the produced pions 
satisfy, approximately, the (generalized) Bj or ken- 



Gottfried condition [24 . I25I ]. whereby the four- 
momentum of a produced particle and the space- 
time position at which it is produced are linearly 
related. Such a correlation between space-time 
and momentum-energy is also a feature of the 
Lund string model, which, incorporated in Jet- 
SET [14], is very successful in describing detailed 
features of the hadronic final states of e + e~ an- 
nihilation. 

A model which predicts such a Q- 
dependence while incorporating the Bjorken- 
Gottfried condition is the so-called r-model, in- 
troduced in ref. [2a]. In this model, it is as- 
sumed that the average production point in the 
overall center-of-mass system, x = (t,r x ,r y ,r z ), 
of particles with a given four-momentum p = 
(E,p x ,p y ,p z ) is 



x»{p^ 



In the case of two-jet events, a 
m t 



yjm 2 +pl 



VE 2 



arpf 1 . (4) 

1/mt where 
2 is the transverse 



P 



mass and r = W t* — r 2 is the longitudinal proper 
time. For isotropically distributed particle pro- 
duction, the transverse mass is replaced by the 
mass in the definition of a and r is the proper 
time. In the case of three-jet events the relation is 
more complicated. The second assumption of the 
r-model is that the distribution of x^(p^) about 
its average, 5a(x^(p^) — x^(p^)), is narrower than 
the proper-time distribution, H(t). The emission 
function of the r-model is 



S(x,p) 



&tH(t)5&(x - arp)pi{p) , (5) 



where H(t) is the (longitudinal) proper-time dis- 
tribution, the factor 5a (x — arp) describes the 
strength of the correlations between coordinate 
space and momentum space variables and p\(p) 
is the experimentally measurable single-particle 
spectrum. 

The two-pion distribution, P2(pi,P2), is re- 
lated to S(x,p), in the plane- wave approximation, 
by the Yano-Koonin formula [23]. The resulting 
two-particle Bose-Einstein correlation function is 
indeed found to depend on the invariant relative 
mometum variable Q, as well as on the values of 
a of the two particles [28]: 



R 2 {pi,P2) = 1 + ReH 



H 



a 2 Q 2 



(6) 
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where H{u) = J 6ltH(t) exp(iur) is the Fourier 
transform (characteristic function) oiH(r). Note 
that H(t) is normalized to unity. This formula 
simplifies further if R2 is measured with the re- 
striction a\ ~ a% ~ a. In that case, for two-jet 
events R2 becomes 



%</>,./'_>)-• 1 +AR(>/7- ( J^- 



(7) 



i?2 is found to depend only on the invariant rela- 
tive momentum Q. Further, the model predicts a 
specific dependence on a, which for two-jet events 
is a specific dependence on Wit [3^ ] . 

Since there is no particle production be- 
fore the onset of the collision, H(t) should be 
a one-sided distribution. We choose a one-sided 
Levy distribution [2^] , which has the characteris- 
tic function j2^] (for a/1) 



Thus for a given average of a of the two particles, 



H{uj) = exp 



(At|o;|) ( 1 — i sign(w) tan 



/air 
V~2~ 



+ i OJTq 



(8) 



r 



where the parameter tq is the proper time of the 
onset of particle production and At is a mea- 
sure of the width of the proper-time distribution. 
For the special case a = 1, see, for example, [30]. 



We have tested that parameter tq is within errors 
zero, hence we fixed it to zero. Using this sim- 
plification, the characteristic function in Eq. © 
yields 



i? 2 (Q, ai, a 2 ) = 7 < 1 + A cos 



tan 



(?) 



r 



exp 



ArQ 5 



1? + Oo 



The result of fitting Eq. © to the l3 two-jet 
event Bose-Einstein correlation data is presented 
on Fig.[TJ The best values of the fit parameters 
and their errors are shown in Table HI Their first 
uncertainty is statistical, the second systematic. 
The confidence levels are shown in Table [TT] for 
varying the transverse mass of the particles inde- 
pendently. The results indicate, that the r- model 
is consistent with the l3 data, the fit quality is 
good and the model is able to describe data well 
not only the low relative momentum region, but 
also the region of 0.5 < Q < 1.6 GeV, where 
anti-correlations are observed. In the large rela- 
tive momentum region of Q > 1.6 GeV, no sig- 
nificant long-range correlations are found and the 
corresponding parameter e is measured to be zero 



within errors. Based on the analysis of the Bose- 
Einstein correlations and the single particle spec- 
tra in terms of the r-model the space-time evolu- 
tion of the particle emitting source can also be re- 
constructed [28]. The first L3 preliminary results 
on such an extremely fast movie were reported in 
refs. [20|, l2l. 
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TABLE I: Results of the fit of Eq. © for two-jet 
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systematic. 
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